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The fluid equations for the baryon-electron system in an expanding universe are derived 
from the Boltzmann equation. The effect of the Compton interaction is taken into account 
properly in order to evaluate the photon-electron collisional term. As an application, the 
acoustic motions of the baryon-electron system after recombination are investigated. The 
effective adiabatic index 7 is computed for sound waves of various wavelengths, assuming 
£N| ' the perturbation amplitude is small. The oscillations are found to be dumped when 7 

C^l \ changes from between 1 (for an isothermal process) to 5/3 (for an adiabatic process). 
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, One of the biggest interest in astrophysics is to understand how structures in the universe are formed. 

0^ ' Observations reveal that the universe is not simply homogeneous and isotropic but contains rich structures 

from stars to the large scale structure such as clusters of galaxies, super-clusters, the great wall and voids 
in it. Although the gravity itself does not have any particular scale, several physical scales associated 
with the structure formation, such as the Jeans scale, the Silk damping scale and the horizon scale at 
the matter-radiation equality epoch, are derived by the inclusion of an effect of primeval radiation |^ |3| . 
Among them, the Jeans scale determines whether small density fluctuations can grow against pressure 
or not. 

The linear evolution of small density perturbations has been well understood. Outside the Jeans scale, 
their evolution is described by the growing mode solution . Assuming the spherical symmetry, we can 
follow the evolution of over-dense regions. The critical threshold of a density contrast to collapse is 1.69 
in the Einstein-de Sitter universe and we can estimate the fraction of the collapsed mass by employing 
the extrapolation of the linear perturbation theory Q . If the wave length of fluctuations is smaller than 
the Jeans scale, however, fluctuations cannot grow up to the critical threshold but begin to oscillate as an 
acoustic wave. Therefore an investigation of the time evolution of the Jeans scale is crucial in particular 
for the study what is the first collapsed object in the universe. 

The Jeans scale is ruled by the sound speed c s and the free fall time of the object or the expansion time 
of the universe. After the recombination of hydrogen atoms, the radiation pressure becomes ineffective 
and equations of state are described by the adiabatic index 7 of baryonic matter. If the energy transfer 
between the baryon and the photon is efficient, we expect isothermal process, i.e., 7 = 1. However, if the 
time scale of the energy transfer becomes longer than the oscillation time scale of the sound wave, the 
adiabatic 7 = 5/3 must be achieved. These two 7's give about factor two difference for the Jeans scales 
in mass. This point is crucial when we investigate the evolution of the baryon density perturbations 
smaller than the Jeans scale after recombination. Moreover, the temperature of gas is very close to the 
radiation temperature at the recombination epoch but gradually separates from it. Eventually, it does 
change the dependence on the red-shift z from 1 + z to (1 + z) 2 . In order to obtain an accurate Jeans 
scale, we need to take into account all these details. In this paper, we formulate the energy transfer 
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between baryons, electrons and photons from the first principle, i.e., the Boltzmann equation. Similar 
treatment on photons has been done by Hu, Scott & Silk and Dodelson & Jubas ||. These works 
focused on the perturbations of photons which are massless particles in order to get anisotropics of the 
cosmic microwave background radiation. Here we formulate perturbation equations of baryon-electron 
fluid which is the massive particle system. The energy transfer between this fluid system and photons is 
coming through the Compton scattering. We study in detail the sound speed and the Jeans scale after 
the recombination. It is found that damping of the sound wave is efficient during the change of 7. 

In §11, we describe the Boltzmann equation in the perturbed expanding universe. In §111, the equations 
of the baryon-electron system are derived based on the fluid approximation, by integrating the Boltzmann 
equation for the momentum space. The integrations of the collisional term which describes the Compton 
interaction between electrons and photons is summarized in §IV. The perturbation equations of the 
baryon-electron fluid are obtained in §V. In §VI, we investigate the acoustic motion of the baryon- 
electron fluid after the recombination in the expanding universe as an application of our formalism. §VTI 
is devoted to summary and discussions. In appendix A, we consider the perturbation equation for the rate 
equation in order to complete the perturbed equations obtained in §V. In appendix B, we put summary 
of the physical scales for the cosmological baryon perturbations in the universe. The evolution of the 
matter temperature is summarized in appendix C. 

We will work in units where c = Ti — ks = 1. 



II. BOLTZMANN EQUATION 

We write the perturbed space-time to the Friedmann-Robertson- Walker space-time in the Newtonian 
gauge as 

\ 2 

a 



ds z = g^dx^dx" = -(1 + 2*)dr + (_) (1 + 2^)5 ij dx l dx> , (2.1) 

with introducing the perturbed gravitational potential \& and the curvature perturbation <I>. a is the scale 
factor and suffix indicates the present value. Note that we employ a flat orthogonal coordinate system 
besides the scale factor as indicated by the Kronecker's delta Sij for a background space. \& satisfies the 
Poisson equation 

\ 2 

a 



= 4.Gp(-j S , (2.2) 

where p and 5 are total background density and density fluctuation, respectively, and $ = — "J when the 
anisotropic stress is negligible, e.g., in the matter-dominated era. 

We write next the Boltzmann equation for the distribution function ft a \ (t, x l , q l ) of the (a)-particle as 

df {a) df (a) dx l df (a) dq l _ 

~dT + ^~dt + - c[S{a)] ' (2 - 3) 

with the collision term in the right hand side. This Boltzmann equation is written in terms of the 
momentum q % measured by an observer in the cosmological rest frame. In order to solve this equation, we 
must rewrite the terms dx l / dt and dq l / dt in terms of x l and q % , which is obtained from the equation of 
motion of the (a)-particle. To obtain this equation, we write the 4-momentum in the locally orthonormal 
frame as (q°, q l ). Here the energy in this frame q° is defined as 



<•(<*) 

where q 2 = J^iil 1 ) 2 - 



g°:=,/q 2 + m (a) 2 , (2.4) 
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The equation of motion is given from the geodesic equation. However, the geo desic equation of the 
(a)-particle is commonly written in terms of the 4-momentum p^ in the frame ( |2.l[ ), which is defined by 
p^ = dx^/dX. Here X(d\ = ds/mi a \) is the afhne parameter. We use the fact that the 4-momentum 
(q°, q 1 ) is related to as follows, 



dx 1 a 

dX ao 

?° = v/=^^ = (l + *)p° ■ 



Equations ( |2.5| ) and (2J3) give the following relations, up to the first order of ^ and <E>, 



dx 1 
~~dt 



5r = ?a+*-*)i. 



(2.5) 
(2.6) 

(2.7) 



and 



dg ,; 



<9$ <9$ oV 



ao V dt dxi dt ) 



x j \ a 
IT IP + - 



ao at 



where the over-dot denotes f-differentiation. 

On the other hand, the geodesic equation in the leading order of the perturbation derives 



(2.8) 



dp 1 
~dt 



a a 



a 
ao 



q>q % 



(1 - $)q* + 2 — <S>q l + 2$,j— + 5 l] ( $ J? ° - <i> 



-.2 \ 1 



(2.9) 



ao a ao q" 

where t i = d/dx 1 . Hereafter we will omit 5^ and will write as e.g., ^ = S^^j. Inserting this to Eq.( |2.8| ). 



dq i 
~dt 



ao 
a 



a a a a u g u 



(2.10) 



Now we can write down the left hand side of the Boltzmann equation ( |2.3| ) using the equations fl2.7| ) and 
( [2.10 ). If we employ the conformal time defined by (a/ao)dr) = dt instead of the proper time, it becomes 
as 



^k + ^k ( i_$ + vi/)^ 

dr/ dx 1 q° 



df(a) 

' dq l 



-q l - $Y - ^^q° - j 



2 1 



ao 



c '[/(«)]> 



(2.11) 



where the prime denotes ^-differentiation. Until now, we have assumed only the smallness of perturbations 
to the homogeneous background. 

If we introduc e a f urther assumption that the motion of the (a)-particle is non-relativistic, the Boltz- 
mann equation ( |2.1l| ) reduces to 



df(a) d/(a) <j_ dfja) 

drj dx 1 mioA dq l 



<&')q l — $jm 



(«) 



go 



C[/(a)] 



(2.12) 



where v l = q 1 /rn/ a ^ and terms of 0(v 2 x or <&) and 0(v 3 ) are omitted. This is a familiar non-relativistic 
Boltzmann equation in the expanding universe, when the term proportional to $' is neglected. 
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III. FLUID APPROXIMATION 



We consider a system of non-relativistic particles (baryons and electrons) and photons which interact 
only with the electrons. Neutral and ionized hydrogen atoms and neutral helium atoms are taken into 
account as baryonic components. We further take a single fluid approximation for this baryon and electron 
system since the time scale of the interaction between them is s hort enough. Under these assumptions, 
we here derive the fluid equation from the Boltzmann equation (2.12). 

Since the particles are non-relativistic, we take the distribution function of the (a)-particle, where 
a = e, H and He for the electron, the hydrogen and the helium, as 



/(a) = Sll( a )(x) 

with a normalization of 



2tt 



m (Q )T b (x) 
d 3 q (Q) 



3/2 



exp 



-(q(q) - TO (a) V b (x)) 2 

2TO (Q )T b (x) 



f( a ) = Sn( Q )(x) 



(3.1) 



(3.2) 



Here ?b(x) and v b (x) are the temperature and the peculiar velocity of this fluid system which the suffix 
b denotes. Then it follows that 



d 3 q( Q ) . . 
(2tt) 3 g O)-fr a) = m (c) Sn ( a )Vh, 

^2^)3 ] q U)1( a )f(a) = m { a ) 2 Sn {a) v h l v h 3 + m {a) Sn {a) T h 5 VJ , 



^ 2 ^3 } q(a) 2 g( a )/(a) = m (a) 3 Sn {a) v h 2 v h l + 5to (q) 2 Sn {a) T h v h 



(3.3) 
(3.4) 
(3.5) 



Operating the follo wing integration and the summation for the Boltzmann equation (2.12) with the 
distribution function (3.1), 



(a) 



(2tt) 3 



we obtain the continuity equation 

dpb 
drj 

for 



d ( 

d^{ PhVh 



Pb(x) = y^m( a )Sn( g )(x). 

(a) 

The collision terms should cancel out by the summation. 



(3.6) 



(3.7) 



(3.8) 



Operating the following integration and the summation to Eq.(2.3) 



E 

(a) 



d 3 q (a 



(«) 



(2tt) s 



(3.9) 



we get the Euler equation 
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where 



and 



___ + 4^- + $ j Pb , b + — ^ + — + Ph — _ AV Compton , (3.10) 



P = J2Sn (a) T h , (3.11) 

(a) 



AT/ 1 . . 

v Compton — 

a 



a /' <i 3 q( 



/ 27r ^3 9(c) C [/(c)] Compton • ( 3 -12) 



The contributions from the collision term between the baryon and the electron should cancel out while the 
contribution from the Compton interaction between the electron and the photon, AVcompton, remains. 
The explicit form of AV2 ompton is given in the next section. 



Operating the following integration and the summation to Eq.(2.3) 



(a) 

we get the energy equation 



dP (a! w \ „ 9* 

d-q 



where 



and 



- + $'JP + Q^Vb % pb = AS Compt o„ , (3.14) 

5P _ 5j2( a ) Sn (a) T b 



2 Pb 2 J2{a) m (a)Sn( a) 



(3.15) 



AE'compton = J (2^P ^^[/(e)] Compton ■ (3.16) 



IV. COLLISION TERM 



Here let us evaluate the collision term between the electron and the photon. The explicit form of 

C[f(e)} Compton IS 



C[/(c)(q)]Compton — 

1 f f f 2d 3 p 2d 3 p' 2d 3 q' 



■{2n) 4 6 {4) {p + q-p' -q')\M\ 



2q° J J J (2ir) 3 2p° (2ir) 3 2p'° (27r) 3 2g'° 
x|(l + / (7 )(p))/ (7 )(p')/(c)(q') - (1 + / (7 )(p'))/( 7 )(p)/(c ) (q)}, (4.1) 

where /( 7 )(p) is the photon distribution function, p and p' are the photon momenta, p° = |p|, p'° = |p'| 
and 5^\p + q - p' - q') = 5(p° + q° - p'° - q'°)S < - 3 ^(p + q p' q') is the 4-dimensional Dirac's delta 
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function. \M\ is the matrix element summing and averaging over the electron spin and the photon 
polarization [H] which is described as, 



M 



(47T) 2 2 
o "EM 



P_ P_ 

P P' 



sin 2 /? 



(4.2) 



Here p and [3 are the photon energy and the photon scattering angle, respectively, in the electron rest 
frame, and a^ M = 3m e 2 OT /Sn. 

Since the explicit integration to get the collision term is very tedious, instead of performing the inte- 
gration of (4.1), we here consider the Boltzmann equation for the photon distribution function 



d/( 7 ) , d/( 7 ) dx i df h) dp i 
dt dx % dt dp 1 dt 



= C[/( 7) ] 



Compton 



(4.3) 



where 



C[/ (7) (P)] Compt> 

l 

2p° 
1 



(2n) 4 5 [A) {p + q-p' -q')\M\ 



2d 3 q 2d 3 q' 2d 3 p' 
(27r) 3 2g° (27r) 3 2</ (2ir) 3 2p'° 



X-<j (1 + /( 7 )(p))/( 7 )(p')/(c)(q , ) - (1 + /( 7 )(p'))/(7)(P)/(c)(q) }■ 



(4.4) 



Note the factor of 1/2 in th e las t line. We should put this factor because we have defined the electron 
distribution function by Eq.(3.1) or Eq.(4.5) , wh ere the spin degree of freedom is summed. 

The integration of the collision term Eq.( [4.4|) has been performed by Dodelson & Jubas || and Hu, 
Scott & Silk Rj. The calculation is based on the assumption that the electron motion is non-relativistic, 
and is performed by expanding in terms of 0(l/m ). In particular, in the paper by Dodelson & Jubas 
B , the integrat ion of the collision term is carried by assuming the same form of the electron distribution 
function as Eq.(3.1), i.e., 



/ (e ) = n e (x) 



2tt 



m e Tb(x) 



3/2 



exp 



(q c - m v b (x)) 



2 1 



2m e T b (x) 



(4.5) 



and by expressing the photon distribution function in the power expansion of 0(l/m e ) up to the second 
order 



/( 7 ) = /( 7 ) (0) b ) + /( 7 ) (1) (p) + /( 7) (2) (p)- 

Following their result || R the collision term is written in the second order of perturbation as 

^/(7)(P)]compton = C«+C( 2 \ 

with 



n c a T 



/(7)| ) 1) + ^(7)^ ) P 2 (M)-/(7) (1) -P ^^ 



(4.6) 



(4.7) 



(4.8) 



and 



*Eq.(5.1) in their paper seems to contain a typographical error. 



G 



c« = 



n c a T 



/(7)^+^/(7)fP 2 (M)-/( 7 ) (2) 



(0) 



■w^^-o* 8 + 1) + «b V) 2 ^h^m 2 + 



(<9p° 



3 

20' 



,(0)' 



where 



/(7)F 



i 2 



Vb ■ P 

l V b||p| 



with P;(/i) being a Legendre polynomial. 
The energy transfer rate is described as 



Ai?Compton — 2 

a 



J^^P° C \-hl) (P)]compto 



Therefore we obtain 

A-Ecompton 



a 4n e CTTAy(o) / to, 



a rn 



^v b ^(x)-T 7 (°'+T b (x) 



(4.9) 

(4.10) 
(4.11) 



(4.12) 



where T^ ^ is the background photon temperature and p 7 
energy density. 

The momentum transfer rate can be found from 



(0) 



(4.13) 

ir 2 T 7 ^ 4 /l5 is the background photon 



= 2 



d 3 p 



PC[/(7)(p)]compto 



/ Compton — " I , n \q 

a a J (2ir) 6 

and the result of the integration in the leading order of perturbations gives the well known form 



AVcompton = \ — rtcOT/4 0) 

3 a 



(v 7 - v b ) . 



(4.14) 



(4.15) 



V. PERTURBATION EQUATIONS 

Elcctrons(e), neutral and ionized hydrogen atoms(H) and helium atoms(He) are the particle species of 
the baryon-electron fluid. The number densities for each species are written 

n e = x(l - )n B , (5.1) 



2 

?i H = f 1 - Vpjn-B, (5.2) 

"He = j"b, (5.3) 

where y p is the primordial helium mass fraction, njj and njjo are the number density of neutral and 
ionized hydrogen and helium, respectively, ub = «h + 4riH C is the total baryon number density and x is 
an ionization fraction defined as x = n c / (njj + 2nn e )- Then, 
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Pb = y^,m( a )Sn( a ) ~ m p n B 

(a) 



p=J2 Sn («) T * = 1 

(a) V 



■ a; • 



y P 



(5.4) 
(5.5) 



where m p is the proton mass. 

First let us summar ize the equations for the baryo n-electron fluid in the expanding universe. The 
continuity equation is (3.7). The Euler equation ( 3 . 1 0| ) reduces to 



+ — + $' )v b i +v h 



;dvh 1 dP 3* a n e cTT 
dxi pb dx 1 dx 1 ao R 



(5.6) 



with the use of Eq.(3.7) and R = 3p\y/4pj. The energy equation (3.14) reduces to 



dr,\ 2 

dx' 1 



■Vb 



•(i-y P /2 



Vb 



(0) 

a T p\ 



Vb 



d f Vb 

dx 1 



h - 



1 dP _5_/a 
Pb di] Pb^a 



(7 + *)' 



T(°)-r b (x)-^vb ; 



by using Eqs.(|3.7|), (|5.l|) and (5.4). We supplement the equation of state, 



_5_P 

2p h 



P = 



Ph 



l + x 



Vv 



b • 



(5.7) 



(5i 



Now we solve the above equations by a perturbative method assuming a small deviations from the 
uniformity. Define the perturbative expansions as follows, 



Pb = Pb ( %)(l + £b(?7,x)), 
h = hW( v )(l + A h ( V: x)), 

P = P(°>(7 7 )(l + Ap(7 y ,x)), 

7b = Tb (0) W(l + AT b (77,x)), 
x = x ^ (77) + fa(77,x), 



together with Vb = Vb(?y,x) 



The zero-th order equations of (3.7) and (5.7) are 



Pb 



(o)' + 3 « (o) = o 



(5.9) 
(5.10) 
(5.11) 
(5.12) 
(5.13) 



(5.14) 



and 



1 (dP<® 



+ 2— 



5~P<°> = 4 



,.(()) 



(l - lfe/2 



(o) 
OT/O7 



(Tf -T b (0) ), (5.15) 



respectively. Eq.(5.15) reduces to 1 1 1 



(o)' + 2 ^ Tb (o) 



i (r (o)_ Tb ( 0))j 



(5.16) 
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where we have neglected the term proportional to x^'T^ ^ [j] and used the zero-th order equation of 
state, 



ft<°> = 



5 p(o) 



P (o) = e^l( 1 + x m _ yp Uo) + s\ \ Tb ( 0) _ 

rn p \ 2 V 2 / / 



2p b (°) ' 

We defined the Compton energy transfer time scale rys as 

3a z(°)(l-y p /2)a x p^ 



(0) 



3 a 



mjl + 5c(°) - (ar(°) + 3/2)y p /2\ 



The solution of Eq.( 5.16 ) is discussed in appendix C. 
The perturbed equations of (3.7) and (5.6) are 



5 h ' + 3$' 



dvb 
dx l 



0. 



and 



t'b + — Ub* + — THT^r— + tt- = 5- K 

a pb 1 - J COT ox 1 a R \ ' 



Vb 



The perturbed equation of (5.7) yields 



5 2 p(°)' 

3 ft(°) ~ 3 P(°) 



r(0) 



Tb (0) 



1 + 1(0) - (as(0) + 3/2)y p /2 
where we have used the perturbative part of the equation of state as follows, 

A h =A P -S h , 

and 



A P = 6 b + A a 



1 + - (a;(°) + 3/2)y p /2 



(5-17) 



(5.18) 



(5.19) 



(5.20) 



(5.21) 



(5.22) 



(5.23) 



In order to complete the perturbation equations, w e nee d an e quat ion to specify the evolution of the 
perturbation of the ionization rate Sx in equations ( 5.21 ) and ( 5.23 ), which we have considered in the 
appendix A. In the early phase of the recombination the terms proportional to Sx become important. 



' As we will see in the below, the matter temperature follows the photon temperature in the recombination 
regime. This is because the energy transfer time scale tje is small enough. After recombination, the fraction of 
the residual ionization x w is almost fixed, and the time variation of x<V is small. This is the reason why we 
neglected the term proportional to x^°''T^ '. 
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VI. SOUND WAVES AFTER THE RECOMBINATION 



As an application of the basic relations obtained in the previous sections, we consider the behavior of 
the sound wave of the baryon-electron fluid after the recombi natio n epoc h. Negl ecting a perturbation of 
the gravitational potential, we get basic equations from Eqs.( [5.19 ), ( 5.20 ), ( 5.21 ) as 



r / , dv h l 



Vb H Vb 

a 



P(°) dA F 



= 0, 



and 



A' 



5 h<-°y 2 p(°y 

3 hi°) 3 P(°) 



;5 h ' = -% x A fc . 



(6.1) 
(6.2) 

(6.3) 



Here we have also neglected a perturbation of the ionization fraction and an interaction with the photon 
in the Euler equation because the momentum transfer rate is so small after the decoupling The wave 
equation is derived from Eqs.(B.l) and (6.2), together with Eq.(5.22), as 



a' P (0) d 2 

5b" + ~*b' - — ToT 7T2 (Ah + *b) = 0. 



(6.4) 



This must be coupled with the energy equation (x3) in order to investigate the sound oscillation of the 
baryon-electron fluid. 

Now let us omit the time-dependence of the background quantities. When the ti me s cale o f th e sound 
oscillation is shorter than the Hubble time, this is a good approximation. Then Eqs.(6.4) and (6.3) reduce 
to simple equations 



8b(k, rf)" + -^k 2 (A h (k, rf) + 5 h (k, rj)) = 0, 
A h (k, rf)' = ^S h (k, v )' - ^A h {k, 7]), 



(6.5) 
(6.6) 



where we took the Fourier mode expansion by setting 5b — db(k,r])e l]5L ' x , and Ah — Ah(k, 77)e* k ' x . The 
above equations yield 



'/E 



d f d 2 5b{k,ri) 
drj \ drj 2 



c f k 5b{k,rf) 



d 2 (5 b (fc,7?) 
dr/ 2 



c k Sb{k,rj) 







Here we defined 



5 P(°> 



Cf 



3 pb 



(o)' 



p(0) 

— (or 



(6.7) 



(6.8) 



As is well known, Cf is the sound speed for the adiabatic state and c c is the one for the isothermal state. 
Taking the wave solution Sb(k,rj) oc e~ luJV , we get the following dispersion relation, 

- ir?Ew(cj 2 - cf 2 k 2 ) +uj 2 - c 2 k 2 = . (6.9) 

In order to solve this equation, it is convenient to introduce the variable w such as ui = iw and we have 

2 5 s 



+ 3 + 



ZD \ 2 

c c k/ 



1 = 0. 



(6.10) 
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The sound speed c s and the adiabatic index 7 are defined as 

l = — 2 - (6.11) 

Therefore we need to solve the cubic equation (|6.10[ ) in order to get the sound speed of the baryon-electron 
system. 



Tm 



FIG. 1. The behaviors of 7 as a function of 1/(1 + 2:). Each lines for the mass scales M = 1M , 1O 3 M , 10 6 M Q , 
1O 9 M0, 1O 12 M0, respectively, with the cosmological parameters, h — 0.5, Qo = 1-0, 51b = 0.1. 7 is changed in 
the earlier stage from 1 to 5/3 for the smaller scale perturbations. 

A cubic equation can be solved exactly in an analytic form while the expression is complicated. It may 



be instructive to show the solution of Eq.(B.lO) in a simple form with some approximation. Expanding the 
solution around the adiabatic state, i.e., expanding in terms of e by setting the solution oj/k — Cf(l + e), 
we find 

gggfl- *±ffg* V (6.12) 

k \ 5(1 + (kctljE) 2 ) J V ' 

It is interesting to clarify the behavior of the adiabatic index 7 in the expanding universe after the 
recombination, which we demonstrate as an example of the usefulness of our formalism in the below. As 



is clear from Eq.(6.9), the adiabatic index 7 is ruled by the ratio of the Compton energy transfer time 
scale ?7e to the sound oscillation time scale. If t]euj <C 1, u> ~ c e k and 7 ~ 1. On the other hand, if 
i]E^ S> 1, u> ~ Cfk and 7 ~ 5/3. Therefore the sound speed depends on scales of the perturbations, and 
it is expected that 7 is changed from 1 to 5/3 as the universe expands. 



Fig.l shows 7 as a function of 1/(1 + z) by directly solving Eq.(3.1C). The ionization fraction is 



calculated properly by solving the recombination process in the expanding universe with the cosmological 
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parameters in the figure caption |Tl],[T^]. As is expected, 7 is changed from 1 to 5/3 in the earlier stage 
of the universe for smaller scale perturbations. 



FIG. 2. The physical mass scales for the baryon perturbations. The cosmological parameters taken here are 
same as those in Fig.l The definitions of these scales are summarized in appendix B. In the figure, 'Horizon' 
means the horizon scale, 'Silk' does the diffusion damping scale, 'tight couple' does the breaking scale of the tight 
coupling approximation of baryon and photon fluids, 'Jeans' does the Jeans scale. 



To get a physical insight, we refer to the familiar illustrative figure, Fig. 2, which gives temporal vari- 
ations of various physical sizes in the expanding universe. The definitions of the curves are summarized 
in appendix B. (See also ref. p3| .) The adiabatic index 7 after the recombination is changed when the 
Compton energy transfer time is equal to the sound oscillation time scale. In this figure, we show the line 
on which the two time scales are equal, i.e., r/^kcc = 1. Note that we use the (baryon) mass scale in the 
unit of the solar mass instead of k, with employing the relation M — (47r/9b/3)(7ra//cao) 3 . The mass and 
the comoving wave number which satisfy the relation rj-Ekc c = 1 are expressed (see also appendix B), 

9.2x^(1 - yj2) , . o 1 , 

fc = rm + Mpc" 1 , 6.13 

( M 3 T W )1/2 

4.6 x 10 10 Oi 3 T( 0) ) 3 / 2 , x q 2 , , 

= 2(0)3(1-^/2)3 (1 + Z) ^ M0 < (6 ' 14) 

where we set the temperature of the microwave background at present T^°\to) — 2.726K, Tb is the baryon 
temperature in unit of Kelvin, and fx = 1 — 3y p /4. 

As is shown in Fig. 2, the Jeans scale after the recombination has the plateau. In this stage the energy 
transfer between background photons and the baryon fluid is effective through the residual electrons, and 
the matter temperature follows the photon temperature. As the universe expands, however, the energy 
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transfer time rises above the Hubble expansion time. After that epoch the matter temperature cools 
adiabatically and drops as Tb oc 1/a 2 . The broken corner of the plateau is the critical time that the 
two time scales become equal, where this epoch is roughly estimated as z ~ 1000(ft b /i 2 ) 2 / 5 [0 (see also 
appendix C). 



FIG. 3. The behaviors of fx> as the function of 1/(1 + 2). Each lines are for the mass scales, M = IMq, 10 3 Mq, 
10 6 Mq, 1O 9 M0, 1O 12 M0, respectively. The cosmological parameters are same as those in Fig.l The damping 
phenomenon occurs when the 7 is changed from 1 to 5/3. 



The line of r]Ekc c = 1 crosses at the broken corner of the plateau of the Jeans mass scale. This 
necessarily happens because of the following reason. The Jeans scale is the scale at which the sound 
oscillation time is equal to the free fall time of the perturbation or the Hubble expansion time. Since the 
broken corner of the plateau is the epoch when the Hubble expansion time is equal to the energy transfer 
time through the Compton interaction between the background photons and the baryon-electron fluid, 
the cross point of two lines is the point when the sound oscillation time, the Hubble expansion time and 
the energy transfer time become all same. 

Now let us discuss an interesting result derived from the dispersion relation ( |6.9| ). As it is apparent 
from the approximated solution, eq.( 3.12 ), u> generally has an imaginary part. This imaginary part of u 
represents the exponential damping of the wave oscillation if Im to < 0. Since the period of the sound 
oscillation is T = 2n/\ Rew|, the damping factor for the amplitude during one period is 



f D = exp 



-2tt 



Im lu 



Re u 



(6.15) 



We show the behavior of Jd in Fig. 3, in the similar way to Fig.l as the function of 1/(1 + z). As is shown 
in Fig. 3, this damping phenomenon may have an effect on the evolution of baryon perturbations. We 
should be notice that our discussions here are based on the assumption that the solutions have the wave 
form. In other word, this damping process is effective only inside the Jeans scale. 
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FIG. 4. The trajectory of the solution on (Ap— (5b)-plane. The horizontal axis is the baryon density perturbation 
— 5b(k,T]), and the vertical axis is the pressure perturbation Ap(fc, rf). We take r\^kc B — 0.1, 1, 10 and 100. 



In order to understand this damping mechanism, we show another aspect of the solutions of Eq.(6.7). 
Introducing lu-r — Rew, lo\ = Imw, we take the solution 

5 h (k, i)=- cos(ct> R 77) e- Wl ". (6.16) 



The initial value is Si,(k,ri = 0) = —1. Here lu is obtained by solving the dispersion relation Eq.^ 
For r]Ekc c 3> 1 and T]^kc c <C 1, solutions are mere harmonic oscillations with no d ampin g. Therefore the 
solutions around r]Ekc c ~ 1 should be investigated. By employing Eqs. (6.5) and ( [5.22| ), the trajectories 



of the solutions on (Ap — £b)-plane are shown in Fig. 4. The horizontal axis is — <5b(fc, rj) and the vertical 
axis is Ap(k,rj). The non-dimensional quantity rj^kc e is the unique physical parameter of the equation. 
We have chosen r]Ekc c = 0.1, 1, 10, and 100 in Fig. 4. 

From the relation —5b — 6V/V, where V is the volume for unit particle number, the horizontal axis 
can be regarded as the change of volume per unit particle number. Thus we can regard the trajectories 
in Fig. 4 in the similar way to the thermodynamical cycles in (pressure-volume)-plane. The damping of 
the oscillation is in proportion to the deviation in one cycle. The most significant damping occurs when 
r]Ekc e = 1 (Fig. 4(b)). The damping does not happen if the deviation in one cycle is negligible as we can 
see in Fig. 4(d). 
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VII. SUMMARY AND DISCUSSIONS 



We have formulated equations for a fluid system with the electron, neutral and ionized hydrogen atoms 
and neutral helium atoms with taking into account the energy transfer between the background photons 
and the residual electrons through the Compton interaction. Using this formulation, we have studied 
the time evolution of the sound speed and the Jeans scale after the recombination. We found that the 
behavior of the adiabatic index 7 after the recombination is controlled by the ratio of the Compton energy 
transfer time scale to the sound oscillation time scale. Then 7 (or the sound speed) depends on scales 
of the perturbations, and is changed from 1 to 5/3 in earlier stage for smaller scale perturbations. This 
formalism enables us to calculate the linear evolution of the very small scale baryon density perturbations 
II- 

We have also discussed the small damping feature of the sound oscillation when 7 changes from 1 to 
5/3. This effect works inside the Jeans scale and seems to be negligible on scales M 10 6 M©. If early 
rcionization occurs, however, the ionization fraction and the Jeans scale increase. And 7 changes from 
5/3 to 1. The damping of the baryon perturbation would be effective on larger scales during the process 
of the reionization of the universe. It has been pointed out that the neutral gas cloud could have an 
instability in a specific photon background [ fl4| . Since we have not consider the specific situation in our 
paper, an unstable mode does not appear in our equations. In the present paper we have also neglected 
perturbations to the ionization fraction. Future work may be required on this point. 

It will be important to quantify the limitation of our formalism. Our basic assumption is the fact that 
the baryon-electron system is treated as tightly coupled single-fluid. This assumption becomes not being 
correct when the collision time scale of neutral interaction rises above the Hubble expansion time. This 
epoch can be estimated in the following way p"5| . The matter temperature can be written as Tt, = 4.5 x 
10 — 3 (f2b/i 2 ) _2 / s (l+£) 2 |K], after the energy transfer through the Compton interaction becomes ineffective 
(see appendix C). Then the mean free time for collision between neutral hydrogen atoms is t c ~ 1/ [n^av), 
where a is the cross section of neutral atoms a ~ irr B [tb is the Bohr radius), and v = y '3Tb / 'ttih- The 
ratio of t c to the Hubble expansion time is therefore expressed as t c H ~ S.l^oh 2 ) 1 ^ 2 (Qbh 2 )^ 4 / 5 (1 + 
z)- 5 / 2 . Eventually we get the red-shift at t c H = 1 as (1 + z) ~ 1.6 (n h 2 ) - 2 (n b h 2 )-°- 32 . This is small 
enough as long as we consider the linear stage of the density perturbations. 

Numerical simulations and semi-analytic calculations of structure formation employ the linear matter 
power spectrum as their initial conditions. This linear matter spectrum is usually calculated without 
taking into account the baryon pressure term after the recombination. It is appropriate, however, only 
if the structures which are larger than the Jeans scale are considered. According to the hierarchical 
clustering scenario, smaller objects are formed earlier than larger ones. We expect the scale of the first 
collapsing object is very close to the Jeans scale. Therefore very accurate estimate of the matter spectrum 
including the baryon pressure term is necessary to understand the early formation of bound objects fl3|| . 
Once the first collapsing objects are formed, they may cool down through formation of hydrogen molecules 
and may fragment into smaller radiating objects like stars and/or quasars |p|,|l6|-^8| . The first radiating 
objects formed Stromgren sphere around them p9| and may eventually ionize all the surrounding gas 
by UV radiation. This reionization process changes the Jeans scale |2Cj]. After that, full numerical or 
semi-analytic calculations including collapsing objects are required. 
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APPENDIX A: RATE EQUATION 

In order to complete the fluid equations we should add the rate equation. For simplicity, we neglect 
the helium fraction, i.e., tib = nn, and neglect the helium recombination process in the rate equation. 
Then, the rate equation is pi 



lh .( av )± nB c(x 2 -(l-x)-^^), (Al) 

\ I x eq / 



where (crv\ is the rate coefficient for recombination to excited states, x cq is the equilibrium ionization 
fraction which is given by the Saha equation 



C cq 2 _ m P ( '"' I ■> \ ^ p -L:).iu-\ /; 
- Xcq Pb \ 27T 



e -13.6cV/T bj (A2) 



and 



_ l + KAn ls fA , 

Here ni s is the number density of hydrogen in the electron ground state, for which we approximate 
"is = ?ib(1 — x), A is the decay rate from the excited state, (3 e = (av)(m e Tb/2n) 3 / 2 e'~ 3A eV / Te and 
K = (a/a)A^/87r with X a being the Lyman alpha photon wave length. 

The Zero-th order equations of the rate equation and the Saha equation are the same forms of Eqs.( |Al| ) 



and (A2) replaced the variables with the zero-th order quantities. 



Next we consider the perturbation of the rate equation. We define 

x cq = ir oq (0) (?7) + Sx cq (r],x), (A4) 



1G 



then 



dSx a 
drj 



= ±(av)n B ^c\( 

a o L V 



1 - * eq (°) 



C 



2 (1- S (°))l ca (°)'\/^ 

, (o) z \ (l-x^)x W(2-x (0)N I 

" (mi \ - V,- 1 - *^ y-^cq V ^ eq y 



+ 2^°) + -^ 



l - W 0) 



Sx 



(l-x oq (o)) 2 



Sx, 



eq 



(A5) 



where 



-K/3 e n B (0) ((l-a; (0) )<5b-<5a;) 



C (l + ifAn B (°)(l-a;(0)))(l + if(A + /3 e )nB( )(l-a;( ))) ' 



We have assumed that (cry) and (3 e are constant. 

The perturbation of the Saha equation can be written as 



dx cq (p(°\ Ph ^) 



0Pb<°> 



P<0) 



phiQ)Sh+ 8x eq{ P^) 



dP(o) 



P (o) A 



p- 



AO) 



(A6) 



(A7) 



APPENDIX B: PHYSICAL SCALES 



In this appendix, we summarize the physical scales which are important for the evolution of the baryon 
perturbations on small scales. The definitions of the physical scales in Fig. 2 are given. Mass scales are 
defined by the amount of baryonic components inside the systems. Here we have set that the temperature 
of the microwave background at present T 7 (i ) = 2.726A'. We write /„ as the neutrino fraction of the 
energy density in the massless particles, and /„ = 0.405 in case of the standard three families of massless 
neutrinos. 

First of all, we use the notation for the physical wave number (length) and the comoving wave number 
(length) as 

fc comv = J^ fcP hy S; A comv = ^^phys = 2 n / k COmV . (Bl) 

It will be useful to give the relation between the red-shift z and the scale factor normalized at the matter 
radiation equality a cq , 

— = 4.04 x 10 4 (1 - U)n h 2 (l + z)- 1 . (B2) 

&eq 



1. Horizon Scale 



We define the horizon wave number and the horizon mass as 

1 



Lcomv =V, (B3) 
H 



phys \ 3 

(B4) 
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which derive 



k£ mv = 3.35 x KT 2 (^1 + ^-- lj (1 - f^^noh 2 Mpc" 1 , 

m h = 9.57 x io 17 (./i + — - 1) (i - u)-^ 2 n h h 2 (n h 2 y 3 M Q . 



2. Jeans scale before the recombination 

We define the Jeans wave length (wave number) and the Jeans mass as 



>phys = 27T = / 7TC S 2 



Mj 

where 



(T) 










U? hys ) 



c 2 ~ 1 



3(1 + i2) ' 

and p m = pb + Pdm with pdm being the energy density of the dark component. Then we have 

2\ 1/2 



kT mv = 1-42 x lCT 1 + (1 - AO^r) (1 + R ) l ' 2 ( l - .M 1/2 Voh 2 Mpc" 1 , 
Mj = 

/ 2N-3/2 

1.25 x W^p + (i - M^r) (i + R)- 3/2 (i .ur 3/2 ^oh 2 r 3 n h h 2 m . 



3. Jeans scale after the recombination 

We can define the Jeans scale after the recombination as 



with 

p(0) 



^phys 



c 2 = 7 



Pb<°> 



= 9.18 x 10- 14 7M7b (0) , 



where /x = 1 — 3y p /4 for :r/ ) <C 1, and Tb is the matter temperature in unit of Kelvin. Then \ 

/ \-i/2 

fc comv = 2 n x 1Q 5 7;uTb _^ ) (! _ f^l/2 Qoh 2 MpQ -^ 

V a eq/ 

(\ 3/2 
7M T b — ) (1 - U)- 3 / 2 {Q h 2 )- 3 Q h h 2 M Q . 
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4. Diffusion damping scale 

We define the diffusion damping scales as |Q 



1 V 1 f 1 j ^ 2 +4(l + i?)/5 



fccomvy 6 J T 1 (1 + i?) 2 



n m2 . (B16) 



An f tt \ 



3 



M D = — Pb — t— , (B17) 



Then we have 



(B18) 



where f = n e <jT(a/cio). For R < 1, we get 

jbg™ = 1-38 x 10 2 ( ^ (0 ff; q) /2) ) - /,) 5/ W 2 ) 3/2 W 2 ) 1/2 Mpc- 1 , 

M D = 1.38 x m^-^0^y\l - yg-iB/^^-Q/a^-i/a Mq) (B19) 

where it(y) = (VTTy(16 - 8y + 6y 2 ) - 16)/15. 

5. Breaking Scale of the tight coupling approximation 

The breaking scale of the tight coupling approximation is defined by ^/k^ s = l/n e <jT, i.e., 

k^T = n e a T -, (B20) 



fcgr = 3.77 x ^( ^(oS^f/a) ) ^ 1 " Z,) 2 W 2 ) 2 (^ 2 ) Mpc" 1 , (B22) 
M BR = 6.75 x io- 1 ^ -^^— ^ (1 - U)-\n a h 2 )- & {n h h 2 r 2 M . (B23) 



6. 7 transition epoch after the recombination 

The adiabatic index 7 after the recombination is changed when the Compton energy transfer time is 
equal to the sound oscillation time scale. We define this scale by l/k^° mv = c c j]e, which leads 

K - (/i 3 Tc) i/ 2(a/acq)3 Mpc , (B24) 

Mc - xW3(l- ^2)3(1 -/„)» M Q , (B25) 

where Tb is the matter temperature in unit of Kelvin as is mentioned above. 
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APPENDIX C: MATTER TEMPERATURE 



As we have derived in §V, the matter temperature follows 

T b (°" + 2-r b (°)^ ??E 1 (rf -r b (°>), 



(Cl) 



where is the photon temperature, and the Compton energy transfer time scale ije is defined as 



t _8fl x(°)(l-y p /2)a T p( ( 



(0) 



3 a o mJ 1 + x(°) - (x(°) + 3/2)y p /2 



The formal solution is 



T b W = 1 jf a 2 Tf r^ 1 exp [- ^ dr/'ry E 1 



(C2) 



(C3) 



We can obtain the epoch when the matter temperature deviates from the photon temperature as follows. 
This epoch is naturally defined as the epoch when the function r]^ 1 exp(— J vo dij'r]^ 1 ^ takes its maximum 

value. Thus we need to solve the equation r]^ 1 + (r]^ 1 ) 2 = 0. From Eq. (|C^), 



Ve 



1 = 1.8 x 10 8 (l-/„) 3 (SVi 2 ) 3 



x^(l-y p /2) ( a 



i - 3y P /4) 



— Mpc- 1 . 



Therefore we get the epoch as 



2.3 x 10 3 x^ 5 (n h 2 ) 4 / 5 , 



(C4) 



(C5) 



' cq 



with neglecting the time variation of the ionization fraction x^°\ Here we set y v = 0.23 and /„ — 0.405. 
We approximate the fraction of the residual electron as pM 



(C6) 



Then from Eq.(C5), we conclude that the matter temperature decouples from the photon's at (1 + z) ~ 
1000(f2 b /i 2 ) 2 / 5 . After this epoch, the matter temperature is adiabatically cooling. 

According to the fully numerical calculation p"2j , the matter temperature in this adiabatic cooling 
phase is well reproduced by the formula, 



= 4.5 x 10- 3 (1 + z) 2 {U h h 2 )- 2 ' b K. 



(C7) 
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The fluid equations for the baryon-electron system in an expanding universe are derived 
from the Boltzmann equation. The effect of the Compton interaction is taken into account 
properly in order to evaluate the photon-electron collisional term. As an application, the 
acoustic motions of the baryon-electron system after recombination are investigated. The 
effective adiabatic index 7 is computed for sound waves of various wavelengths, assuming 
the perturbation amplitude is small. The oscillations are found to be dumped when 7 
changes from between 1 (for an isothermal process) to 5/3 (for an adiabatic process). 
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I. INTRODUCTION 

One of the biggest interest in astrophysics is to understand how structures in the universe are formed. 
Observations reveal that the universe is not simply homogeneous and isotropic but contains rich structures 
from stars to the large scale structure such as clusters of galaxies, super-clusters, the great wall and voids 
in it. Although the gravity itself does not have any particular scale, several physical scales associated 
with the structure formation, such as the Jeans scale, the Silk damping scale and the horizon scale at 
the matter-radiation equality epoch, are derived by the inclusion of an effect of primeval radiation [1-3]. 
Among them, the Jeans scale determines whether small density fluctuations can grow against pressure 
or not. 

The linear evolution of small density perturbations has been well understood. Outside the Jeans scale, 
their evolution is described by the growing mode solution [4,5]. Assuming the spherical symmetry, we can 
follow the evolution of over-dense regions. The critical threshold of a density contrast to collapse is 1.69 
in the Einstein-de Sitter universe and we can estimate the fraction of the collapsed mass by employing 
the extrapolation of the linear perturbation theory [6]. If the wave length of fluctuations is smaller than 
the Jeans scale, however, fluctuations cannot grow up to the critical threshold but begin to oscillate as an 
acoustic wave. Therefore an investigation of the time evolution of the Jeans scale is crucial in particular 
for the study what is the first collapsed object in the universe. 

The Jeans scale is ruled by the sound speed c s and the free fall time of the object or the expansion time 
of the universe. After the recombination of hydrogen atoms, the radiation pressure becomes ineffective 
and equations of state are described by the adiabatic index 7 of baryonic matter. If the energy transfer 
between the baryon and the photon is efficient, we expect isothermal process, i.e., 7 = 1. However, if the 
time scale of the energy transfer becomes longer than the oscillation time scale of the sound wave, the 
adiabatic 7 = 5/3 must be achieved. These two 7's give about factor two difference for the Jeans scales 
in mass. This point is crucial when we investigate the evolution of the baryon density perturbations 
smaller than the Jeans scale after recombination. Moreover, the temperature of gas is very close to the 
radiation temperature at the recombination epoch but gradually separates from it. Eventually, it does 
change the dependence on the red-shift z from 1 + z to (1 + z) 2 . In order to obtain an accurate Jeans 
scale, we need to take into account all these details. In this paper, we formulate the energy transfer 
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between baryons, electrons and photons from the first principle, i.e., the Boltzmann equation. Similar 
treatment on photons has been done by Hu, Scott & Silk [7] and Dodelson & Jubas [8]. These works 
focused on the perturbations of photons which are massless particles in order to get anisotropies of the 
cosmic microwave background radiation. Here we formulate perturbation equations of baryon-electron 
fluid which is the massive particle system. The energy transfer between this fluid system and photons is 
coming through the Compton scattering. We study in detail the sound speed and the Jeans scale after 
the recombination. It is found that damping of the sound wave is efficient during the change of 7. 

In §11, we describe the Boltzmann equation in the perturbed expanding universe. In §111, the equations 
of the baryon-electron system are derived based on the fluid approximation, by integrating the Boltzmann 
equation for the momentum space. The integrations of the collisional term which describes the Compton 
interaction between electrons and photons is summarized in §IV. The perturbation equations of the 
baryon-electron fluid are obtained in §V. In §VI, we investigate the acoustic motion of the baryon- 
electron fluid after the recombination in the expanding universe as an application of our formalism. § VII 
is devoted to summary and discussions. In appendix A, we consider the perturbation equation for the rate 
equation in order to complete the perturbed equations obtained in §V. In appendix B, we put summary 
of the physical scales for the cosmological baryon perturbations in the universe. The evolution of the 
matter temperature is summarized in appendix C. 

We will work in units where c = Ti = ks = 1. 



II. BOLTZMANN EQUATION 

We write the perturbed space-time to the Friedmann-Robertson- Walker space-time in the Newtonian 
gauge as 

ds 2 = g^dx^dx" = -(1 + 2^)dt 2 + (^—^j (1 + 2(S>)8 ij dx i dx j , (2.1) 

with introducing the perturbed gravitational potential \? and the curvature perturbation $. a is the scale 
factor and suffix indicates the present value. Note that we employ a flat orthogonal coordinate system 
besides the scale factor as indicated by the Kronecker's delta 6{j for a background space, \? satisfies the 
Poisson equation 

V 2 * = teGp(^j 8 , (2.2) 

where p and 6 are total background density and density fluctuation, respectively, and $ = — \? when the 
anisotropic stress is negligible, e.g., in the matter-dominated era. 

We write next the Boltzmann equation for the distribution function f( a ){t, x % , q l ) of the (a)-particle as 

df (a) df (a) dx { df (a) dq { _ 

-^r + ^^r + ^^- c[f ^ ] ' (2 - 3) 

with the collision term in the right hand side. This Boltzmann equation is written in terms of the 
momentum q* measured by an observer in the cosmological rest frame. In order to solve this equation, we 
must rewrite the terms dx % j dt and dq % j dt in terms of x % and q % , which is obtained from the equation of 
motion of the (a)-particle. To obtain this equation, we write the 4-momentum in the locally orthonormal 
frame as (q° , q l ). Here the energy in this frame q° is defined as 



q° ■= ^ 2 + m (a f, (2.4) 

where q 2 = Y.M % ) 2 ■ 
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The equation of motion is given from the geodesic equation. However, the geodesic equation of the 
(a)-particle is commonly written in terms of the 4-momentum p^ 1 in the frame (2.1), which is defined by 
p^ = dx^/dX. Here X(dX = ds/m^ a j) is the affine parameter. We use the fact that the 4-momentum 
(q°,q i ) is related to pf , as follows, 



, dx % a ,„ , 

/gTi-rr = - i + $ 

dX do 



q° = V-ffoo 



dx° 
~dX 



Equations (2.5) and (2.6) give the following relations, up to the first order of 'f and $, 



dx l p l 
~dt ~ p° 



(2.5) 
(2.6) 

(2.7) 



dq % a ■ a /9$ 9$ dx 1 \ ■ a . . dp 1 

~JT = 'I + — «T + TTT^T P + — 1 + $ hf ' 
ar a ao \ ar / ao ar 

where the over-dot denotes ^-differentiation. 

On the other hand, the geodesic equation in the leading order of the perturbation derives 



dt 



2--(l - $) g < + 2 -$ g i + 2$ j + & ( * jq° ~<fj\ 

a a a q ' \ q ' 



q 



(2.8) 



(2.9) 



where )8 - = d/dx* . Hereafter we will omit 8' 3 and will write as e.g., 'f 8 - = Inserting this to Eq.(2.8), 

,2 " 



dt 



ao 
a 



a a i a k i * 9 1 

q l + — $g' + 

ao a ao </ 



(2.10) 



Now we can write down the left hand side of the Boltzmann equation (2.3) using the equations (2.7) and 
(2.10). If we employ the conformal time defined by (a/ao)dr) = dt instead of the proper time, it becomes 
as 



dq l 



q 



$y ^-^Af- 



a 



(2.11) 



where the prime denotes ry-differentiation. Until now, we have assumed only the smallness of perturbations 
to the homogeneous background. 

If we introduce a further assumption that the motion of the (a)-particle is non-relativistic, the Boltz- 
mann equation (2.11) reduces to 



df(a) df (a) q i df (a) 
drj dx 1 m( a ) dq l 



a 



C[f {a) ] 



(2.12) 



where v % = </V m (a) an d terms of 0(v 2 x 'f or $) and 0(v 3 ) are omitted. This is a familiar non-relativisti 
Boltzmann equation in the expanding universe, when the term proportional to $' is neglected. 
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III. FLUID APPROXIMATION 



We consider a system of non-relativistic particles (baryons and electrons) and photons which interact 
only with the electrons. Neutral and ionized hydrogen atoms and neutral helium atoms are taken into 
account as baryonic components. We further take a single fluid approximation for this baryon and electron 
system since the time scale of the interaction between them is short enough. Under these assumptions, 
we here derive the fluid equation from the Boltzmann equation (2.12). 

Since the particles are non-relativistic, we take the distribution function of the (a)-particle, where 
a = e, H and He for the electron, the hydrogen and the helium, as 

3/2 r "(q(«) " m (a) v b (x)) 2 ] 



/ W = 5n W (x) (m (a ^ b (x)) 6XP 



2m (a )T b (x) 



with a normalization of 

j3 



J^-f(a) = Sn (a) ( X ) . (3.2) 

Here 7b(x) and v b (x) are the temperature and the peculiar velocity of this fluid system which the suffix 
b denotes. Then it follows that 



^yr9(a)/(«) = m (a) Sn (a) v h \ (3.3) 

9(a) ?(«)/(«) = m {a) 2 Sn {a) v h l v h 3 + m {a) Sn {a) T h S l: i , (3.4) 

"(«) 'lit 3 c 2 i i r 2 o rr i /o c\ 

(2?r) 3 q(a) q (a)J(a) = m(a) i«(a)V b v h + 5m (a ) Sn (a) T h v h . (3.5) 



Operating the following integration and the summation for the Boltzmann equation (2.12) with the 
distribution function (3.1), 

(a) V > 

we obtain the continuity equation 



for 



Pb(x) = ^ra( a )5n (a )(x). (3.8) 

(a) 



The collision terms should cancel out by the summation. 

Operating the following integration and the summation to Eq.(2.3), 



E/ $f «!■>■ (»■»> 



we get the Euler equation 



4 



where 



__ + 4^- + $jp b , b +— ^ + ^-+ Pb ^--AT/ Compton , (3.10) 



P = £Sn (a) T b , (3.11) 

(a) 



AV , = SL /" rf3q ( e ) 

v Compton — 

a 



^2^3 9(e) C '[/(e)]compton • (3.12) 



The contributions from the collision term between the baryon and the electron should cancel out while the 
contribution from the Compton interaction between the electron and the photon, AVcompton, remains. 
The explicit form of AVc ompton is given in the next section. 

Operating the following integration and the summation to Eq.(2.3), 

(a) > 

we get the energy equation 

BP fa' \ 9* 

—^--5 -+$' )P+TT- «bVb = AVcompton , (3-14) 

or] V a / ox 1 



where 



and 



> S f= ^°> S "<f , (3,5, 



A-Ecompton = ~ / TTTTT qf e ^[/(e)]compton • (3.16) 

2m e a J (2TT) d { > y > 



IV. COLLISION TERM 

Here let us evaluate the collision term between the electron and the photon. The explicit form of 

C[/(e)] Compton 1® 



C[/(e)(q)]compton — 

1 f f f 2d 3 p 2d 3 p' 2rf 3 q' 



-(2ir) 4 8( 4 )(p + q-p' -q')\M\ 



2q° J J J (2tt) 3 2p° (2tt) 3 2p'° (27r) 3 2 g '° 

X {(1 + / (7 )(p))/( 7 )(p')/(e)(q') - (1 + /( 7 )(p'))/(7)(P)/(e)(q)}. (4-1) 

where /( 7 )(p) is the photon distribution function, p and p' are the photon momenta, p° = |p|, p'° = |p'| 
and h^ A \p + q — p' — q') = 6(p° + q° — p'° — q'°)6^ 3 \p + q — p' — q') is the 4-dimensional Dirac's delta 
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function. M is the matrix element summing and averaging over the electron spin and the photon 
polarization [9] which is described as, 



M 



(47T) 2 



'EM 



p 



p 



p ft 



■ sin 2 (] 



(4.2) 



Here p and [3 are the photon energy and the photon scattering angle, respectively, in the electron rest 
frame, and aj| M = 3m e 2 <7T/87r. 

Since the explicit integration to get the collision term is very tedious, instead of performing the inte- 
gration of (4.1), we here consider the Boltzmann equation for the photon distribution function 



df(j) 9/( 7 ) dx i 5/ (7 ) 



dt 



dx l dt dp 1 dt 



C[f {l) ] 



7) JCompton 1 



(4.3) 



where 



C[/( 7 )(p)]compt 

V / / J (27r) 3 2g° (27r) 3 2 g '° (27r) 3 2p'° 
1 



2rf 3 q 2rf 3 q' 2rf 3 p' 



-^^(p + q-p' -q')\M\ 



X 2 Ul + /(7)(P))/(7)(p')/(e)(q') " (1 + /( 7 )(p'))/(7)(P)/(e)(q) 



(4.4) 



Note the factor of 1/2 in the last line. We should put this factor because we have defined the electron 
distribution function by Eq.(3.1) or Eq.(4.5), where the spin degree of freedom is summed. 

The integration of the collision term Eq.(4.4) has been performed by Dodelson & Jubas [8] and Hu, 
Scott & Silk [7]. The calculation is based on the assumption that the electron motion is non-relativistic, 
and is performed by expanding in terms of 0(l/m e ). In particular, in the paper by Dodelson & Jubas 
[8], the integration of the collision term is carried by assuming the same form of the electron distribution 
function as Eq.(3.1), i.e., 



/(e) = «e(x) 



27T 



m e T b (x) 



3/2 



exp 



-(q e - m e v b (x)) 2 
2m e Xb(x) 



(4.5) 



and by expressing the photon distribution function in the power expansion of 0(l/m e ) up to the second 
order 



/( 7 )=/( 7 ) (0) (p°) + /( 7 ) (1) (p) + /( 7 ) (2) (p). 
Following their result [8] *, the collision term is written in the second order of perturbation as 

C[/( 7 )(p)]compton =C (1 '>+C (2 \ 

with 



n e a T 



/(7)^+ 2^ JP ^)"/(7) 



(D_pO 



5/ (7) 



(o) 



dp 



-/iv h 



(4.6) 



(4.7) 



(4.8) 



*Eq.(5.1) in their paper seems to contain a typographical error. 



n e a T 



/(7)^ ) + ^(7)^P 2 M-/(7) (2) 



W/^^0< 2 + i) + »bV) 2 



dp 



d' 2 f< S 0) 11 3 

(dp ) 2 4(T 20 ; 



5/( 7) (0) 



- 1 d/j 



(oy 



where 



ff ^ 



WW , 



Vb P 
| v b||p| 



with P;(/u) being a Legendre polynomial. 
The energy transfer rate is described as 



AEi 



Therefore we obtain 



AE, 



Compton — ^ I / ^ \ q 

a J (2TT) d 
a 4n e aTP y (o) ( m, 



P C[/( 7 )(p)] 

Compto 



Compton 



0>o nif 



(4.9) 

(4.10) 
(4.11) 



(4.12) 



where T 7 ( -°- ) is the background photon temperature and pj 
energy density. 

The momentum transfer rate can be found from 



^^v b 2 (x)-T 7 (°)+T b (x)), (4.13) 
7r 2 T 7 ( -°- )4 /15 is the background photon 



AVc, 



P 



I'Compton — ~~ L J ^2~^p"PC[/(7)(p)]compton , 

and the result of the integration in the leading order of perturbations gives the well known form [5], 

(4.15) 



(4.14) 



4 a 



AVcompton = 77— "eO-T/cl - 1 I V 7 - V b 

3 a y \ ' 



V. PERTURBATION EQUATIONS 

Electrons(e), neutral and ionized hydrogen atoms(H) and helium atoms(He) are the particle species of 
the baryon-electron fluid. The number densities for each species are written 

n e = x(l - y)«B, (5.1) 
(l-y P )n B , (5.2) 

"He = ^~ n B, (5-3) 

where y p is the primordial helium mass fraction, «h and njje are the number density of neutral and 
ionized hydrogen and helium, respectively, «b = n H + 4«He is the total baryon number density and x is 
an ionization fraction defined as x = n e / («h + 2«He)- Then, 
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Pb = ^2 m (a) Sn (a) - m p n B (5.4) 

(a) 

P = ^5n (a )T b = (l + x-^(x + ^)jn B T h , (5.5) 

(a) V 7 

where m p is the proton mass. 

First let us summarize the equations for the baryon-electron fluid in the expanding universe. The 
continuity equation is (3.7). The Euler equation (3.10) reduces to 

dv h l fa' A • :dv h l 1 dP 9* a n e a T j ■ -\ 
with the use of Eq.(3.7) and R = ?>pb/ Ap^\ The energy equation (3.14) reduces to 



drj \ 2 J \ a / \ 2 / 9a; 8 \ 2 J Pb dr) pb V a 

9* a ^(l-yp/2)«7 T ^ 0) / me X 

+ 3^ b = 4 ^ ^ ^H-T^x)-^ j, (5.7) 

by using Eqs.(3.7), (5.1) and (5.4). We supplement the equation of state, 

"=£( i+ *-?(* + §)K 

Now we solve the above equations by a perturbative method assuming a small deviations from the 
uniformity. Define the perturbative expansions as follows, 

Pb=Pb (0) (r?)(l + ^b(r?,x)), (5.9) 

/ l = / l (°)(^)(l + A ft (r y ,x)), (5.10) 

F = p(°)( I) )(l + Ap( I) ,x)), (5.11) 

T b =T b (°)(r y )(l + A Tb (r y ,x)), (5.12) 

x = x(°\rj) + 6x(r 1 ,x), (5.13) 

together with v b = v b (ry,x). 

The zero-th order equations of (3.7) and (5.7) are 

Pb (0) ' + 34 b (°) = 0, (5.14) 

and 

^ + 2 a -h^ _ + ,°L P (A = 4^ " (0) ( 1 -^ /2 )" T ^°V T (0) _ Tb (oA > (5 . 15) 
dr/ a Pb^ \ dr) a J ao rn e m p V 7 /' 

respectively. Eq.(5.15) reduces to [10] 

Tb (°)' + 2 ^T b (°) = r,^(T^ - T b (°)), (5.16) 
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where we have neglected the term proportional to x^'T^ ^ t and used the zero-th order equation of 
state, 

*<". = ^, P<.» = ^( 1+I .»»-|( I »4))T t <». (5.17) 

We defined the Compton energy transfer time scale rj^ as 

3 "o mjl + x(°) - (x(°) + 3/2)t/ p /2j 

The solution of Eq.(5.16) is discussed in appendix C. 
The perturbed equations of (3.7) and (5.6) are 

V + 3$' + ^=0, (5.19) 



and 



dx l 



,./ a 1 ,. P(°) dA P 5* a n e a T / ,- , , 



The perturbed equation of (5.7) yields 



A/ /5M°)' 2P(°)V 2 £/ 
~*~ ' 3 TiW 3 ~PW" ) h ~ 3 b 



- -La , ** 4 0) -r b <°> , (i-yp/a)^ \ 

" % ^ h + xW Tb ( °) 1 + *(°) - (*(°) + 3/2)y p /2j ' 1 j 

where we have used the perturbative part of the equation of state as follows, 

A h =A P -8 b , (5.22) 

and 

Ap = S b + A Th + 1 + , (0) V _ ( , ( o )+ ; 3/%p/2 • (5-23) 

In order to complete the perturbation equations, we need an equation to specify the evolution of the 
perturbation of the ionization rate 6x in equations (5.21) and (5.23), which we have considered in the 
appendix A. In the early phase of the recombination the terms proportional to 6x become important. 



*As we will see in the below, the matter temperature follows the photon temperature in the recombination 
regime. This is because the energy transfer time scale rj-E is small enough. After recombination, the fraction of 
the residual ionization x^ - 1 is almost fixed, and the time variation of x^ - 1 is small. This is the reason why we 
neglected the term proportional to x^'T^°K 
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VI. SOUND WAVES AFTER THE RECOMBINATION 



As an application of the basic relations obtained in the previous sections, we consider the behavior of 
the sound wave of the baryon-electron fluid after the recombination epoch. Neglecting a perturbation of 
the gravitational potential, we get basic equations from Eqs.(5.19), (5.20), (5.21) as 

V + ^ = 0, (6.1) 



a p b 



/5 hW' 2 P(°)'\ 2 
^+(1^5-3^)^-3^ = -^. < 6 - 3 ) 

Here we have also neglected a perturbation of the ionization fraction and an interaction with the photon 
in the Euler equation because the momentum transfer rate is so small after the decoupling The wave 
equation is derived from Eqs.(6.1) and (6.2), together with Eq.(5.22), as 

/ p(0) f)2 

6 b " + -6 b ' - —--^(A h +6 b ) = 0. (6.4) 

This must be coupled with the energy equation (6.3) in order to investigate the sound oscillation of the 
baryon-electron fluid. 

Now let us omit the time-dependence of the background quantities. When the time scale of the sound 
oscillation is shorter than the Hubble time, this is a good approximation. Then Eqs.(6.4) and (6.3) reduce 
to simple equations 

«b( W + —^k 2 (A h (k,Ti) + 6 b (k, n)) = 0, (6.5) 

A h (k, rj)' = h b (k, V )' - V ^A h (k, V ), (6.6) 

where we took the Fourier mode expansion by setting 6 b = 6 b (k , ry)e 8k * x , and = Ah(k, ry)e 8k * x . The 
above equations yield 

'4( ! T + c ' V4t,M) ) + = ° • < 6 ^ 

Here we defined 

2 5P(°) 2 

Cf =3^)' Ce =^M- (6 ' 8) 

As is well known, Cf is the sound speed for the adiabatic state and c e is the one for the isothermal state. 
Taking the wave solution 8 b (k,rj) oc e~ Ujjrl , we get the following dispersion relation, 

-iry E w(w 2 - Cf 2 k 2 ) +lo 2 - c 2 k 2 = . (6.9) 

In order to solve this equation, it is convenient to introduce the variable w such as u> = izu and we have 



c e krjE 



(5)((5) 2+ 5) + (5) 2+ — (610 » 
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The sound speed c s and the adiabatic index 7 are defined as 



a, = — Im 



7 



(6.11) 



Therefore we need to solve the cubic equation (6.10) in order to get the sound speed of the baryon-electron 
system. 




FIG. 1. The behaviors of 7 as a function of 1/(1 + 2). Each lines for the mass scales M 
10 9 Mq, 10 12 Mq, respectively, with the cosmological parameters, h = 0.5, fio = 1.0, 
the earlier stage from 1 to 5/3 for the smaller scale perturbations. 



1M , 1O 3 M , 1O 6 M , 
= 0.1. 7 is changed in 



A cubic equation can be solved exactly in an analytic form while the expression is complicated. It may 
be instructive to show the solution of Eq.(6.10) in a simple form with some approximation. Expanding the 
solution around the adiabatic state, i.e., expanding in terms of e by setting the solution u>/k = c f (l + e), 
we find 



ui 
k 



c f [ 1 



1 + ikc f r] E 
5(1 + {kc fm f) 



(6.12) 



It is interesting to clarify the behavior of the adiabatic index 7 in the expanding universe after the 
recombination, which we demonstrate as an example of the usefulness of our formalism in the below. As 
is clear from Eq.(6.9), the adiabatic index 7 is ruled by the ratio of the Compton energy transfer time 
scale r/E to the sound oscillation time scale. If tjeoj < 1, w ~ c e k and 7 ~ 1. On the other hand, if 
rjE^J ^> 1, w ~ Cfk and 7 ~ 5/3. Therefore the sound speed depends on scales of the perturbations, and 
it is expected that 7 is changed from 1 to 5/3 as the universe expands. 

Fig.l shows 7 as a function of 1/(1 + z) by directly solving Eq.(6.10). The ionization fraction is 
calculated properly by solving the recombination process in the expanding universe with the cosmological 
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parameters in the figure caption [11,12]. As is expected, 7 is changed from 1 to 5/3 in the earlier stage 
of the universe for smaller scale perturbations. 
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FIG. 2. The physical mass scales for the baryon perturbations. The cosmological parameters taken here are 
same as those in Fig.l The definitions of these scales are summarized in appendix B. In the figure, 'Horizon' 
means the horizon scale, 'Silk' does the diffusion damping scale, 'tight couple' does the breaking scale of the tight 
coupling approximation of baryon and photon fluids, 'leans' does the leans scale. 

To get a physical insight, we refer to the familiar illustrative figure, Fig. 2, which gives temporal vari- 
ations of various physical sizes in the expanding universe. The definitions of the curves are summarized 
in appendix B. (See also ref. [13].) The adiabatic index 7 after the recombination is changed when the 
Compton energy transfer time is equal to the sound oscillation time scale. In this figure, we show the line 
on which the two time scales are equal, i.e., rj^kc e = 1. Note that we use the (baryon) mass scale in the 
unit of the solar mass instead of k, with employing the relation M = (47T/9b/3)(7ra/A;ao) 3 . The mass and 
the comoving wave number which satisfy the relation rj^kc e = 1 are expressed (see also appendix B), 



9.2a;( )(l - t/p/2)^ 1 
k = \ yp/ ' 1 + z 3 Mpc" 1 , 

(M 0) ) 112 

4.6 x lO 10 ^ 3 ^) 3 / 2 , , „ , 

M= ,(0)3(1-^/2)3 (l + ')- 9 ^ 2 M0, 



(6.13) 
(6.14) 



where we set the temperature of the microwave background at present Tj°\to) = 2.726K, Tb is the baryon 
temperature in unit of Kelvin, and \i = 1 — 3j/ p /4. 

As is shown in Fig. 2, the Jeans scale after the recombination has the plateau. In this stage the energy 
transfer between background photons and the baryon fluid is effective through the residual electrons, and 
the matter temperature follows the photon temperature. As the universe expands, however, the energy 
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transfer time rises above the Hubble expansion time. After that epoch the matter temperature cools 
adiabatically and drops as Tb oc I /a 2 . The broken corner of the plateau is the critical time that the 
two time scales become equal, where this epoch is roughly estimated as z ~ lOOC^f^/i 2 ) 2 / 5 [10] (see also 
appendix C). 
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FIG. 3. The behaviors of fn as the function of 1/(1 + z). Each lines are for the mass scales, M = IMq, 10 3 Mq, 
10 6 Mq, 10 9 Mq, 10 12 Mq, respectively. The cosmological parameters are same as those in Fig.l The damping 
phenomenon occurs when the 7 is changed from 1 to 5/3. 



The line of rjEkc e = 1 crosses at the broken corner of the plateau of the Jeans mass scale. This 
necessarily happens because of the following reason. The Jeans scale is the scale at which the sound 
oscillation time is equal to the free fall time of the perturbation or the Hubble expansion time. Since the 
broken corner of the plateau is the epoch when the Hubble expansion time is equal to the energy transfer 
time through the Compton interaction between the background photons and the baryon-electron fluid, 
the cross point of two lines is the point when the sound oscillation time, the Hubble expansion time and 
the energy transfer time become all same. 

Now let us discuss an interesting result derived from the dispersion relation (6.9). As it is apparent 
from the approximated solution, eq.(6.12), u> generally has an imaginary part. This imaginary part of w 
represents the exponential damping of the wave oscillation if Im u> < 0. Since the period of the sound 
oscillation is T = 2ir/\ Rew|, the damping factor for the amplitude during one period is 



f D = exp 



-2tt 



Im 



Re ui 



(6.15) 



We show the behavior of fo in Fig. 3, in the similar way to Fig.l as the function of 1/(1 + z). As is shown 
in Fig. 3, this damping phenomenon may have an effect on the evolution of baryon perturbations. We 
should be notice that our discussions here are based on the assumption that the solutions have the wave 
form. In other word, this damping process is effective only inside the Jeans scale. 



13 



0. 



1 - 



1 - 



x 



1 - 



J II II | 1 [ II II II 1 II 1 | L 

(a)77 E kc e =0.1 : 


_i i i i i i i i i | i i i i i i i i I. 

I (b)77 E kc e =l : 


_ 1 II II II 1 1 , 

: : 

- i 1 i i i i 1 i i i i i i i i 1 i i i i 1 r 


J II 1 1 II II 1 _ 

- i i i i i 1 i i i i 1 i i i i i i i i 1 r 



-1 -0.5 0.5 1 -1 -0.5 0.5 



FIG. 4. The trajectory of the solution on (Ap- <5b)-plane. The horizontal axis is the baryon density perturbation 
— 8b(k,7]), and the vertical axis is the pressure perturbation Ap(fc,jj). We take r]s,kc e = 0.1,1,10 and 100. 

In order to understand this damping mechanism, we show another aspect of the solutions of Eq.(6.7). 
Introducing = Rew, u>\ = Imw, we take the solution 



&b(k,r)) = - cos(w R ry) e Wir > . 



(6.16) 



The initial value is 8yj(k,r) = 0) = —1. Here u> is obtained by solving the dispersion relation Eq.(6.9). 
For r]Ekc e 1 and rjEkc e <C 1, solutions are mere harmonic oscillations with no damping. Therefore the 
solutions around rj^kc e ~ 1 should be investigated. By employing Eqs. (6.5) and (5.22), the trajectories 
of the solutions on (Ap — (*>b)-plane are shown in Fig. 4. The horizontal axis is —6t,(k,r)) and the vertical 
axis is Ap(k,rj). The non-dimensional quantity rjEkc e is the unique physical parameter of the equation. 
We have chosen rjEkc e = 0.1, 1, 10, and 100 in Fig. 4. 

From the relation —8^ = 6V/V, where V is the volume for unit particle number, the horizontal axis 
can be regarded as the change of volume per unit particle number. Thus we can regard the trajectories 
in Fig. 4 in the similar way to the thermo dynamical cycles in (pressure-volume)-plane. The damping of 
the oscillation is in proportion to the deviation in one cycle. The most significant damping occurs when 
ry E fcc e = 1 (Fig. 4(b)). The damping does not happen if the deviation in one cycle is negligible as we can 
see in Fig. 4(d). 
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VII. SUMMARY AND DISCUSSIONS 



We have formulated equations for a fluid system with the electron, neutral and ionized hydrogen atoms 
and neutral helium atoms with taking into account the energy transfer between the background photons 
and the residual electrons through the Compton interaction. Using this formulation, we have studied 
the time evolution of the sound speed and the Jeans scale after the recombination. We found that the 
behavior of the adiabatic index 7 after the recombination is controlled by the ratio of the Compton energy 
transfer time scale to the sound oscillation time scale. Then 7 (or the sound speed) depends on scales 
of the perturbations, and is changed from 1 to 5/3 in earlier stage for smaller scale perturbations. This 
formalism enables us to calculate the linear evolution of the very small scale baryon density perturbations 

We have also discussed the small damping feature of the sound oscillation when 7 changes from 1 to 
5/3. This effect works inside the Jeans scale and seems to be negligible on scales M 1O 6 M . If early 
reionization occurs, however, the ionization fraction and the Jeans scale increase. And 7 changes from 
5/3 to 1. The damping of the baryon perturbation would be effective on larger scales during the process 
of the reionization of the universe. It has been pointed out that the neutral gas cloud could have an 
instability in a specific photon background [14]. Since we have not consider the specific situation in our 
paper, an unstable mode does not appear in our equations. In the present paper we have also neglected 
perturbations to the ionization fraction. Future work may be required on this point. 

It will be important to quantify the limitation of our formalism. Our basic assumption is the fact that 
the baryon-electron system is treated as tightly coupled single-fluid. This assumption becomes not being 
correct when the collision time scale of neutral interaction rises above the Hubble expansion time. This 
epoch can be estimated in the following way [15]. The matter temperature can be written as Tb = 4.5 x 
10- 3 (fi b /j 2 )~ 2/5 (l + 2) 2 [K], after the energy transfer through the Compton interaction becomes ineffective 
(see appendix C). Then the mean free time for collision between neutral hydrogen atoms is t c ~ l/{n-Q<Tv), 
where a is the cross section of neutral atoms a ~ irrg (re is the Bohr radius), and v = y/3Tb/mn. The 
ratio of t c to the Hubble expansion time is therefore expressed as t c H ~ 3.1(£lo/i 2 ) 1/ ' 2 (^i/i 2 ) _4/ ' 5 (l + 
z)" 5 / 2 . Eventually we get the red-shift at t c H = 1 as (1 + z) ~ 1.6 {£l h 2 )° 2 {Q h h 2 )-° 32 . This is small 
enough as long as we consider the linear stage of the density perturbations. 

Numerical simulations and semi-analytic calculations of structure formation employ the linear matter 
power spectrum as their initial conditions. This linear matter spectrum is usually calculated without 
taking into account the baryon pressure term after the recombination. It is appropriate, however, only 
if the structures which are larger than the Jeans scale are considered. According to the hierarchical 
clustering scenario, smaller objects are formed earlier than larger ones. We expect the scale of the first 
collapsing object is very close to the Jeans scale. Therefore very accurate estimate of the matter spectrum 
including the baryon pressure term is necessary to understand the early formation of bound objects [13]. 
Once the first collapsing objects are formed, they may cool down through formation of hydrogen molecules 
and may fragment into smaller radiating objects like stars and/or quasars [3,16-18]. The first radiating 
objects formed Stromgren sphere around them [19] and may eventually ionize all the surrounding gas 
by UV radiation. This reionization process changes the Jeans scale [20]. After that, full numerical or 
semi-analytic calculations including collapsing objects are required. 
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APPENDIX A: RATE EQUATION 

In order to complete the fluid equations we should add the rate equation. For simplicity, we neglect 
the helium fraction, i.e., «b = "h, and neglect the helium recombination process in the rate equation. 
Then, the rate equation is [10] 

^ = -(o-v)-n B c(x 2 - (1 - x)-^-) , (Al) 
Of] v ' a V l-x eq J 

where (cm) is the rate coefficient for recombination to excited states, x eq is the equiliblium ionization 
fraction which is given by the Saha equation 



X eq 2 _ / m eTb \ ^ g -13.6eV/T b /^2) 

1 - Xeq Ph\2-Kj 

and 

C - 1 + A ' Awis (A3) 

G -l + A-(A + /3 e )n ls • (Ad) 

Here n ls is the number density of hydrogen in the electron ground state, for which we approximate 
n ls = «b(1 — x), A is the decay rate from the excited state, (3 e = (<7i>)(m e Tb/27r) 3/ ' 2 e _3 ' 4 eV / Te and 
K = (a/a)\^/SiT with \ a being the Lyman alpha photon wave length. 

The Zero-th order equations of the rate equation and the Saha equation are the same forms of Eqs.(Al) 
and (A2) replaced the variables with the zero-th order quantities. 

Next we consider the perturbation of the rate equation. We define 

x eq = Xeq^irj) + 8x eq (r],x), (A4) 
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then 



drj do 



where 



[av)n B ^C 



J0) 2 _ MV^\ fsc 



1 - * eq (°) 



c 



+ 6t 



+ 2xW + x °* _ ) 8x - 1 , ' eq )1„ eq > 8x, 



1 " * eq (°) 



(1 " W )) 5 



eq 



8C 



-A'/3 e n B (0) ((l-* (0) )^b-^*) 



C (l + KAn B (°)(l " * (0) )) (l + K ( A + Pe)n B (°Kl ~ * (0) )) 



(A5) 



(A6) 



We have assumed that (cm) and j3 e are constant. 

The perturbation of the Saha equation can be written as 



8x P 



dx eq (p( o ),p b (°y) 



P(o) 



(0) dx eq (pw, Pb (°y) 

Pb 6b + dW) — 



k(°) 



(A7) 



APPENDIX B: PHYSICAL SCALES 



In this appendix, we summarize the physical scales which are important for the evolution of the baryon 
perturbations on small scales. The definitions of the physical scales in Fig. 2 are given. Mass scales are 
defined by the amount of baryonic components inside the systems. Here we have set that the temperature 
of the microwave background at present T 7 (to) = 2.726A'. We write /„ as the neutrino fraction of the 
energy density in the massless particles, and /„ = 0.405 in case of the standard three families of massless 
neutrinos. 

First of all, we use the notation for the physical wave number (length) and the comoving wave number 
(length) as 

fc comv = (— )fc phys , A comv = (— )A phys = 27r/F omv . (Bl) 

It will be useful to give the relation between the red-shift z and the scale factor normalized at the matter 
radiation equality a eq , 

— = 4.04 x 10 4 (1 - f„)tt h 2 (l + z)- 1 . (B2) 



1. Horizon Scale 



We define the horizon wave number and the horizon mass as 

V, (B3) 



Lcomv 
H 



^h = ^P(^) =^P(7i^) . (B4) 



H 
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which derive 



fc™ mv = 3.35x 10 ~ 2 (y^ 1 + ~ " *) t 1 ~ fv) 1/2 V h 2 Mpc" 1 , 
M H = 9.57x 10 17 (^Jl + JL-lj (1- f v )- 3 l 2 tt h h 2 (tt h 2 )- 3 M Q . 

2. Jeans scale before the recombination 

We define the Jeans wave length (wave number) and the Jeans mass as 



,phys _ 2?T _ / 7TC S 2 



where 



c 2 - 1 



3(1 + R)' 

and p m = pb + Pdm with /9<jm being the energy density of the dark component. Then we have 

2\ 1/2 



/ 2 \ i/ z 

= 1-42 x lO" 1 ^ + (1 - /,)^f (1 + /?) 1/2 (1 - Z,) 1 / 2 ^/* 2 Mpc" 1 , 
\ a cr / 

Mj = 

, 2 \-3/2 

1.25 x i i6 ^ + (i _ f v )^-) (1 + i?)" 3 / 2 (l - /,)- 3 / 2 (fi / l 2 )- 3 «b/ l 2 M G 
\ a cr / 



3. Jeans scale after the recombination 

We can define the Jeans scale after the recombination as 



^phys 27T / 7TC, 



fcj phys V °p* 



with 

,(°) 



c * 2 = = 9 - 18 x 10" 14 7Atfb (0) 



where p = 1 — 3j/ p /4 for a;(°) <C 1, and Tb is the matter temperature in unit of Kelvin. Then 

-1/2 



A ,comv = 2 71 x 10 5 / 7AiTb JL) (l-/,) 1 / 2 ^ 2 Mpc" 1 , 

( \ 3/ ' 2 

Mj = 1.81 x 10" 3 7/U T b — (1 - Z,)- 3 / 2 ^/* 2 )- 3 ^/* 2 M . 

V a eq / 
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4. Diffusion damping scale 



We define the diffusion damping scales as [12] 

1 V 1 f 1 , i? 2 +4(l + i?)/5 



fccomvy 6 J f (1 + i?) 2 



/i / \ 3 

47T / 7T 



where f = n e aT(a/ao). For i? <C 1, we get 

-1/2 



*fST = 3-77 x 10 4 ( J*^^ ) " (1 " /,) W) W) Mpc" 1 , 



M BR = 6.75 x 10" 1 (^(o|^2)) (1 " f,r 6 ^oh 2 )- 6 (n h h 2 )- 2 M Q . (B23) 



(B16) 



Md = tHc^' (B17) 



(B18) 



*g>™ = 1.38 x lO ^co"^^) )" (1 " f*) 5/4 (n h*r'\n b h*)W Mpc" 1 , 
M D = 1.38 x 10 7 ( g( or ( ( r-^/2) ) 3/2(1 " /^)" 15/4 («o/ l 2 )- 9/2 ("b/ l 2 )- 1/2 M 0J (B19) 



where = (^/l + t/(16 — 8j/ + 6j/ ) — 16)/15. 

5. Breaking Scale of the tight coupling approximation 

The breaking scale of the tight coupling approximation is defined by l/k^[ s = l/n e CT, i.e., 

kgg v = n e <T T — , (B20) 

do 

M BK=fpb(-i) 3 . (B21) 

^ BR y 

Then we have 

^2 \ -i 



(B22) 



6. 7 transition epoch after the recombination 

The adiabatic index j after the recombination is changed when the Compton energy transfer time is 
equal to the sound oscillation time scale. We define this scale by l/fc™ mv = c e rj^, which leads 

= .^xW /,,,!-/,). 

(// d Te) 1 /^(a/a e q) d 

_ 1.60xl0- 31 (^e) 3 / 2 (a/a eq ) 9 o 2 , 2 ._ 9 , . 

where Tb is the matter temperature in unit of Kelvin as is mentioned above. 
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APPENDIX C: MATTER TEMPERATURE 



As we have derived in §V, the matter temperature follows 

T b (°)' + 2^T b ( ) = ^ 1 (^-T b ( )), (CI) 
where T-y " 1 is the photon temperature, and the Compton energy transfer time scale t/e is defined as 

* = / W ('-^W"' . (CS) 

J "o m e M + x(°) - (x(°) + 3/2)t/ p /2j 



The formal solution is 

1 /•>) r H 

(C3) 



a 



„ . 7 r) E exp 

o 



We can obtain the epoch when the matter temperature deviates from the photon temperature as follows. 
This epoch is naturally defined as the epoch when the function ry^ 1 exp(— dr/ry^ 1 ) takes its maximum 

value. Thus we need to solve the equation jj^ 1 + (jfe 1 ) 2 = 0- F rom Eq. (C2), 

1 = 1.8 x 10 8 (1 - /,) 3 (Q fe 2 ) 3 " (0)(l :^ /2) (-?-) " Mpc" 1 . (C4) 

(l-3t/ p /4) \a eq J 



Therefore we get the epoch as 



a 

a, 



2.3 x 10V°) 2 / 5 (fi o /> 2 ) 4/5 , (C5) 



eq 



with neglecting the time variation of the ionization fraction x^°\ Here we set y p = 0.23 and /„ = 0.405. 
We approximate the fraction of the residual electron as [21] 

*(°) ~ io- 5 (n /» 2 ) 1/2 (nb/» 2 ) _1 . (C6) 

Then from Eq.(C5), we conclude that the matter temperature decouples from the photon's at (1 + z) ~ 
1000(S!b'i 2 ) 2/ ' 5 - After this epoch, the matter temperature is adiabatically cooling. 

According to the fully numerical calculation [12], the matter temperature in this adiabatic cooling 
phase is well reproduced by the formula, 

= 4.5 x 10" 3 (1 + z) 2 (Q b /i 2 )" 2 / 5 K. (C7) 
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